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Abstract 

A simple model with a novel type of dynamics is introduced in order to inves- 
tigate the emergence of self-ordered motion in systems of particles with biologically 
motivated interaction. In our model particles are driven with a constant absolute 
velocity and at each time step assume the average direction of motion of the particles 
in their neighborhood with some random perturbation [r]) added. We present numer- 
ical evidence that this model results in a kinetic phase transition from no transport 
(zero average velocity, |va| = 0) to finite net transport through spontaneous symmetry 
breaking of the rotational symmetry. The transition is continuous since |va| is found 
to scale as {rjc — with (3 ~ 0.45. 



One of the most interesting aspects of many particle systems is that they exhibit 
a complex cooperative behavior during phase transition [1]. This remarkable feature 
of equilibrium systems has been studied in great detail for the last couple of decades 
leading to a deeper understanding of processes which may take place in an assembly 



of interacting particles. Concepts like scaling, universality and renormalization have 
resulted in a systematic picture of a wide range of systems in physics [1,2]. 

Recently there has been an increasing interest in the rich behaviour of systems 
which are far from equilibrium. Processes such as aggregation, viscous flows or biolog- 
ical pattern formation have been shown to involve scaling of the related geometrical 
and dynamic quantities characterizing these phenomena [3,4]. As a further similarity 
with equilibrium systems, the existence of phase transition type behavior has also 
been demonstrated in several investigations of growth processes [5-8] . These analogies 
with the basic features of equilibrium systems have represented a particularly impor- 
tant contribution to the understanding of the complex behavior of non-equilibrium 
processes. 

In this work we introduce a model with a novel type of dynamics in order to in- 
vestigate clustering, transport and phase transition in non-equilibrium systems where 
the velocity of the particles is determined by a simple rule and random fluctuations. 
The only rule of the model is: at each time step a given particle driven with a con- 
stant absolute velocity assumes the average direction of motion of the particles in its 
neighborhood of radius r with some random perturbation added. We show using simu- 
lations that in spite of its simplicity, this model results in a rich, realistic dynamics, 
including a kinetic phase transition from no transport to finite net transport through 
spontaneous symmetry breaking of the rotational symmetry. 

In this sense our model is a transport related, non- equilibrium analog of the fer- 
romagnetic type of models, with the important difference that it is inherently dynamic: 
the elementary event is the motion of a particle between two time steps. Thus, the 
analogy can be formulated as follows: the rule corresponding to the ferromagnetic 
interaction tending to align the spins in the same direction in the case of equilibrium 
models is replaced by the rule of aligning the direction of motion of particles in our 
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model of cooperative motion. The level of random perturbations we apply are in 
analogy with the temperature. 

Beyond the above aspects (analogies and simplicity) the proposed model is 
interesting because of possible applications in a wide range of biological systems in- 
volving clustering and migration. Biological subjects have the tendency to move as 
other subjects do in their neighborhood [9]. In addition to such trivial examples as 
schools of fish, herds of quadrupeds or flocks of flying birds, our model can be applied 
to the less known phenomena during bacterial colony growth [10]. There are bacteria 
(e.g., a strain of Bacillus Subtilis) which exhibit cooperative motion in order to survive 
under unfavorable conditions. The present model with some modifications is already 
capable to reproduce the main observed features of the motion (collective roatation 
and flocking) of bacteria [10]. Other biologically motivated, recent theoretical inves- 
tigations of clustering, aggregation and orientational order in systems with diffusing 
directed objects have concentrated on the possible spatial patterns arising from an 
integro-differential equation approach and from cellular automata type models [11]. 

Furthermore, we expect that our model can be used to interpret the results 
of experiments on clustering and convection in a system of disks floating on an air 
table [12]. These experiments represent a physically motivated possible application 
of the present model, since they are being carried out in order to understand the 
flow of granular materials under specific conditions. We are aware that two groups 
are working on developing models similar to ours in order to interpret these air table 
experiments [13]. 

The actual simulations were carried out in a square shaped cell of linear size L 
with periodic boundary conditions. The particles were represented by points moving 
continuously (off-lattice) on the plane. We used the interaction radius r as unit to 
measure distances (r = 1) while the time unit. At = 1 was the time interval between 
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two updating of the directions/positions. In most of our simulations we used the 
simplest initial conditions: i) at time t = N particles were randomly distributed 
in the cell and ii) had the same absolute velocity v and iii) randomly distributed 
directions 9. The velocities {vj} of the particles were determined simultaneously at 
each time step and the position of the ith particle updated according to 



Here the velocity of a particle Vi(t + 1) was constructed to have an absolute value v and 
a direction given by the angle 9{t + 1). This angle was obtained from the expression 



where {0{t))r denotes the average direction of the velocities of particles (including 
particle i) being within a circle of radius r surrounding the given particle. The av- 
erage direction was given by the angle axctg[{sm{9{t)))r/ {cos{9{t)))r]- In (2) A9 is a 
random number chosen with a unifom probability from the interval [—rj/2, rj/2]. Thus, 
the term represents noise which we shall use as a temperature-like variable. Corre- 
spondingly, there are three free parameters for a given system size: ry, p and v, where 
V is the distance a particle makes between two updatings. 

We have chosen this realization because of its simplicity, however, there may be 
several more interesting alternatives of implementing the main rules of the model. In 
particular, the absolute value of the velocities does not have to be fixed, one can intro- 
duce further kind of interaction between a particles and or consider lattice alternatives 
of the model. In the rest of this paper we shall concentrate on the above described 
simplest version and investigate the nontrivial behavior of the transport properties as 
the two basic parameters of the model, the noise rj and the density p = N/L^ are 
varied. We used v = 0.03 in the simulations we report on because of the following 
reasons. In the limit v ^ the particles do not move and the model becomes an 
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Xi(t + 1) =x,(t) +Vi(t)At. 



(1) 



9{t+l) = {9{t))r + A9, 



(2) 



analog of the well known XY model. For v — > oo the particles become completely- 
mixed between two updates and this limit corresponds to the so called mean field 
behavior of a ferromagnet. We use v = 0.03 for which the particles always interact 
with their actual neighbours and move fast enough to change the configuration after 
a few updates of the directions. According to our simulations in a wide range of the 
velocities (0.003 < v < 0.3) the actual value of v does not affect the results. 

Fig. la-d demonstrates the velocity fields during runs with various selections 
for the value of the parameters p and rj. The actual velocity of a particle is indicated 
by a small arrow, while their trajectory for the last 20 time steps is shown by a 
short continuous curve, (a) At t = the positions and the direction of velocities are 
distributed randomly, (b) For small densities and noise the particles tend to form 
groups moving coherently in random directions, (c) At higher densities and noise the 
particles move randomly with some correlation, (d) perhaps the most interesting case 
is when the density is large and the noise is small; in this case the motion becomes 
ordered on a macroscopic scale and all of the particles tend to move in the same 
spontaneously selected direction. 

This kinetic phase transition is due to the fact that the particles are driven with 
a constant absolute velocity; thus, unlike standard physical systems in our case the 
net momentum of the interacting particles is not conserved during collision. We have 
studied in detail the nature of this transition by determining the absolute value of the 
average normalized velocity 



of the entire system of particles as the noise and the density were changed. This 
velocity is approximately zero if the direction of the motion of the individual particles 
is distributed randomly, while for the coherently moving phase (with ordered direction 
of velocities) Va — l,so that we can consider the average velocity as an order parameter. 




(3) 
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First we gradually decreased the amount of noise rj in cells of various sizes for 
a fixed density p and observed a transition from a disorderly moving phase to a phase 
with a coherent motion of the particles (Fig. 2a). The uncertainity of the data points 
is within the range of the symbols except for runs carried out with 4000 and 10000 
particles close to the transition. For these r] values the statistical errors estimated 
from 5 runs with different initial conditions are in the range of 5% (resulting in an 
overlap of the results for a limited number of rj values) due to the slow convergence 
and large fluctuations. In Fig. 2b we show how Va changes if the noise is kept constant 
and the density is increased. 

Quite remarkably, the behavior of the kinetic order parameter Va is very similar 
to that of the order parameter of some equilibrium systems close to their critical point. 
The strongest indication of a transition in our nonequilibrium model is the fact that as 
we go to larger system sizes the region over which the data show scaling is increasing 
(see Fig. 3a). Only an extremely unusual crossover could change this tendency. A 
plausible physical picture behind our finding is the following: since the particles are 
diffusing, there is mixing in the system resulting in an effective (long range) interaction 
radius. 

Thus, we can assume that in the thermodynamic limit our model exhibits a 
kinetic phase transition analogous to the continuous phase transition in equilibrium 
systems, i.e., 

Va ^ iVcip) - and ~ (p - Pciv))i^ (4) 

where /3 and S are critical exponents and ?7c(p) and pdv) sixe the critical noise and 
density (for L — > oo), respectively. We can determine /? and 5 corresponding to the rate 
of vanishing of the order parameter from plotting logVa as a function of \og[r]c{L) — 
r])/rjc{L)] and log[p — pc{L) / pc{L)] for some fixed values of p and ry, respectively 
(Fig. 3). For finite sizes r]c{L) and Pc{L) are L dependent, thus, we used such values 
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of quantities for which the pfots in Fig. 3 were the straightest in the relevant region 
of noise or density values. The slope of the lines fitted to the data can be associated 
with the critical exponents for which we obtained (3 — 0.45 ± 0.07 and 5 — 0.35 ± 0.06. 
The errors in determining (3 and 5 are due to the uncertainties in the i) Va and the 
ii) r]c{L) and Pc{L) values. Since the scaling plots in Fig. 3 depend sensitively on the 
choice of the critical noise and density and our method of determing their value is 
indirect (from the straightness of the data sets) we give rather conservative estimates 
for the errors of (3 and 5. 

We have carried out a finite size scaling analysis of r]c{L) and obtained ?7c(oo) = 
2.9 ± 0.05 for p = 0.4 (note that the "infinite temperature" limit of our model is 
r/c = 27r). As indicated, r]c depends on p, in fact, we expect a phase diagram (a line 
of critical temperatures) analogous to that of disordered ferromagnets, rjc playing the 
role of temperature and p playing the role of the density of spins. In this case (3 and 5 
are expected to have the same value. On the other hand, strong crossover effects are 
likely to effect their actual values in a finite size simulation. Although our estimates 
for (3 and 5 are different, on the basis of our simulations we cannot exclude the possi- 
bility (allowed by our error bars) that they become equal in the thermodynamic limit. 
However, the determination of the phase diagram and a more precise calculation of 
the exponents of our new model is outside of the scope of the present work which 
concentrates on demonstrating the main features of a novel nonequilibrium system. 

The emergence of cooperative motion in our model has analogies with the ap- 
pearance of spatial order in equilibrium systems. This fact, and the simplicity of 
our model suggests that, with appropriate modifications, the theoretical methods 
for describing critical phenomena may be applicable to the present kind of far-from- 
equilibrium phase transition. The kinetic phase transitions which have been observed 
in surface growth models [5-8] are both in analogy and different from the situation de- 
scribed here. The similarity is in the scaling behaviour of an inherently non-equilibrium 
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order parameter, while the two kinds of processes are distinct from the point of the 
driving force acting on the particles. Self-driven particles are uncommon in physics, 
but they are typical in biological systems, including live organisms and the so called 
"molecular motors" having attracted great interest recently [14]. Transitions have 
been observed in traffic models [15] consisting of particles (cars) which can also be 
interpreted as self-driven particles. 

There are interesting further variations of the model investigated in this work. 
It is expected that taking into account a hard core term in the interaction or using 
semi periodic or open boundary conditions results in additional non-trivial effects. Our 
preliminary results [10] indicate that a model with hard core repulsion and specific 
boundary conditions can be successfully used to interpret recent observations of coher- 
ent motions in geometrically complex bacterial colonies growing on soft agar surfaces 
[10,16-18]. 
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Figure captions 

Figure 1. In this figure the velocities of the particles are displayed for varying values 
of the density and the noise. The actual velocity of a particle is indicated by a small 
arrow, while their trajectory for the last 20 time step is shown by a short continuous 
curve. The number of particles is N = 300 in each case, (a) t = 0, L = 7 77 = 2.0 
(b) For small densities and noise the particles tend to form groups moving coherently 
in random directions, here L = 25, ry = 0.1; (c) After some time at higher densities 
and noise {L = 7 rj = 2.0) the particles move randomly with some correlation; (d) For 
higher density and small noise (L = 5 77 = 0.1) the motion becomes ordered. All of 
our results shown in Figs. 1-3 were obtained from simulations in which v was set to 
be equal to 0.03. 

Figure 2. (a) The absolute value of the average velocity (f^) versus the noise i] in cells 
of various sizes for a fixed density p. The symbols correspond to □ - A?" = 40, L = 3.1; 
+ - N= 100, L = 5; X - N = 400, L = 10; A - N = 4000, L = 31.6. - N= 10000, 
L — 50. In Fig. 2b (for L — 20) we show how Va changes if the noise is kept constant 
and the density is increased. 

Fig. 3 Dependence of logfa on log[(?7c(-^) — 'n)/'nc{L)] and log[(p — pc{L)/ pc{L)]. The 
slope of the lines fitted to the data can be associated with the critical exponents P 
and 5. a) is for p = 0.4, b) is for L = 20 and r] = 2.0 
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